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Abstract. 

Fermi bubbles, the recently observ ed giant (~ 1 kpc high) gamma-ray emit- 
ting lobes on either side of our Galaxy (ISu et al.l l2010). appear morphologically con- 
nected to the Galactic center, and thus offer a chance to test several models of super- 
massive black hole (SMBH) evolutio n, feedback and relation with their host galaxies. 
We use a physical feed back model (King 2003, 20101) and novel numerical techniques 
(Nayaksh in et al.l200 9) to simulate a short burst of activity in Sgr A*, the central SMBH 
of the Milky Way, ~ 6 Myr ago, temporally coincident with a star formation event in 
the central p arsec. We are able to r eproduce the bubble morphology and energetics both 
analytically (Zubovas e t al.ll20llb and numerically (Zubovas & Nayakshin, in prep). 
These results provide strong support to the model, wh ich was also used to simulate 
more extreme environments (iNayakshin & Poweril2010h . 

The AGN radiation radiation pressure drives a wind with a momentum flux_M 0Ut v - 
^Edd/c with v ^ rjc 0.1c, where rj - 0.1 is the radiative efficiency dKingll2003 ). This 



wind shocks against the surrounding gas (perhaps producing y rays) and pushes it away, 
forming an outflow. In the Milky Way, the wind shock cannot cool outside R coo \ ~ 10 pc 
and hence transfers most of the kinetic energy rate (~ 0.05 Z^dd) to the ambient gas (this 
is an energy-driven flo w). Such an outfl ow, while driven, moves with a constant veloc- 



ity v e ~ 1000 km s _1 (|King et al.ll201 ll) . Once the quasar switches off, the shell coasts 
for an order of magnitude longer than the driving phase t q , easily reaching radii of tens 
of kpc. 

The outflow morphology can become non-spherical due to anisotropic matter dis- 
tribution in the Galaxy, such as the dense gas in the Central Molecular Zone (CMZ) 
which is too heavy for even an energy-driven outflow to lift. This qualitatively explains 
the morphology of the Fermi bubbles. We use their observed and inferred properties to 
constrain the gas fraction (ratio of gas density to background potentia l density) in the 



Galaxy halo and the Sgr A* outburst duration (see Zu bovas et al]|201 11 for details) 



/ g <7-10- 3 /; t q > 2.5 • 10 5 yr. (1) 

We test the model numerica lly, using GADGET w ith a 'virtual particle' method 
of implementing wind feedback (Nayaks hin et al.ll2009l) . We embed the SMBH (which 



produces feedback for a time ? q ) and CMZ into a spherically symmetric isothermal halo 
with <x = 100 km/s and a constant f g . We vary the free parameters t q and f g . 

Figure[TJ left, shows that our model, with f q = 1 Myr and f g = 10~ 3 , can repro- 
duce the morphology and size of the observed Fermi bubbles. The CMZ is perturbed 
but not dispersed by the wind and collimates the outflow into two cavities. The total 
energy content inside the cavities is a small fraction of the input and also agrees with 
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Figure 1 . Left: Gas surface density (left) and temperature (right) at t — 6 Myr. The 
outflow is collimated and forms teardrop-shaped cavities with similar morphology to 
that of the observed Fermi bubbles. The simulation has f q = 1 Myr, /_3 = 1. Right: 
Same as left, but f q = 0.3 Myr. The cavity morphology is obviously inconsistent 
with observations. A change in f g produces inconsistent results as well. 



observational constraints. Figure [U right, shows a simulation with t q = 0.3 Myr, which 
produces bubbles clearly inconsistent with observations. We can thus put tight con- 
straints on both parameters. We also require the CMZ mass to be - 10 M n , but its 
aspec t ratio is not important. A physical heating-cooling prescription dSazonov et al. 



2005) does not change the results significantly either. Therefore our findings are quite 
robust with regard to the uncertainties involved in the initial conditions. 

We have shown that our physically motivated SMBH wind feedback model can ex- 
plain the Fermi bubbles. In addition, the same model works for qu asars as their SMBHs 
estab lish the M-cr relation and clear the host galaxies of gas (e.g. [Nayakshin & Power 
2010), suggesting there is no fundamental difference between the processes that were 
responsible for forming the galaxies at z > 2 and the processes what is happening in 
local, mostly dormant, galactic nuclei. 

Acknowledgments. This research used the ALICE High Performance Computing 
Facility at the University of Leicester. KZ is supported by an STFC studentship. 



References 

King, A. 2003, ApJ, 596, L27. |arXiv: astro-ph/Q388342| 

King, A. R. 2010, MNRAS, 402, 1516. Q91 1. 16391 

King, A. R., Zubovas, K, & Power, C. 201 1, MNRAS, L263+. 11(94 . 3682 
Nayakshin, S., Cha, S.-H., & Hobbs, A. 2009, MNRA S, 397, 1314 .[Q9Q5 . 2896 

Nayakshin, S., & Power, C. 2010, MNRAS, 402, 789. (8911 . 24p 

Sazonov, S. Y, Ostriker, J. P., & et al. 2005, MNRAS, 358, 168. arXiv : astr o-ph/Q4 11(886 
Su, M., Slatyer, T. R., & Finkbeiner, D. P. 2010, ApJ, 724, 1044. 1(88 5 .54801 " 
Zubovas, K, King, A. R., & Nayakshin, S. 2011, MNRAS, 415, L21. [1184^5443 



